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Abstract 

A positive cosmological constant simplifies the asymptotics of forever 
expanding cosmological solutions of the Einstein equations. In this pa- 
' per a general mathematical analysis on the level of formal power series 

f"*^ , is carried out for vacuum spacetimes of any dimension and perfect fluid 

. spacetimes with linear equation of state in spacetime dimension four. For 

equations of state stiffer than radiation evidence for development of large 
CO . gradients, analogous to spikes in Gowdy spacetimes, is found. It is shown 

■ that any vacuum solution satisfying minimal asymptotic conditions has a 

""^ | full asymptotic expansion given by the formal series. In four spacetime 

dimensions, and for spatially homogeneous spacetimes of any dimension, 
,* ^ ' these minimal conditions can be derived for appropriate initial data. Us- 

j^Q, ing Fuchsian methods the existence of vacuum spacetimes with the given 

formal asymptotics depending on the maximal number of free functions 
is shown without symmetry assumptions. 



c5 1 1 Introduction 



Spacetimes with accelerated expansion have come to play an important role 
in cosmology. The accelerating phase may be in the early universe (inflation) 
or at the present epoch (quintessence). The simplest way to produce a model 
with accelerated expansion which solves the Einstein equations is to introduce 
a positive cosmological constant. A good survey article on this topic is |T9). 

The fact that a positive cosmological constant leads to solutions of the Ein- 
stein equations with exponential expansion is associated with the term 'cosmic 
no hair theorem'. In the following we investigate possibilities of proving theo- 
rems related to these ideas. In the setting of formal series a satisfactory answer 
is obtained for the Einstein equations in vacuum or in the presence of a perfect 
fluid with linear equation of state. There are formal series solutions which have 
the expected asymptotic behaviour and which depend on the maximum number 
of free functions. This also holds for vacuum spacetimes in higher dimensions. 
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In the case of even space dimensions it is in general necessary to allow terms 
with logarithmic dependence on the expansion parameter. This throws some 
light on what is special about three space dimensions. These results are proved 
in Section 2. 

While most of the results in three space dimensions obtained in Section 2 
confirm the results of one new phenomenon was observed. Evidence is 
obtained that for fluids with an equation of state stiffer than that of a radia- 
tion fluid inhomogeneous structures can be formed. This is reminiscent of the 
formation of spikes near the initial singularity in Gowdy spacetimes |17| . 

In Section 3 it is shown that in the vacuum case minimal assumptions on the 
asymptotics in the expanding phase imply that the spacetime has an asymp- 
totic expansion of the form already exhibited as formal series. Unfortunately 
we do not know in general how to obtain these minimal assumptions starting 
from conditions on initial data. An exception to this is the case of three space 
dimensions where it is shown in Section 4 that the minimal assumptions can 
be deduced from results of Friedrich 0|, [S] based on the conformal method. 
The minimal assumptions can also be verified in the case of certain spatially 
homogeneous vacuum spacetimes of any dimension, as shown in Section 5. In 
particular, there arc genuine solutions of the Einstein equations whose asymp- 
totics contain non- vanishing logarithmic terms. In Section 6 Fuchsian methods 
are applied to show the existence of vacuum spacetimes of any dimension with 
the asymptotics given in Section 2 and depending on the maximum number 
of free functions. Finally, Section 7 shows that in a model problem, the wave 
equation on de Sitter space, full information on asymptotics of solutions with 
arbitrary initial data can be obtained. 

2 Perturbative solutions 

A perturbative treatment of four-dimensional vacuum spacetimes with positive 
cosmological constant can be found in JH]- In that paper formal solutions 
are written down without any mathematical derivation being given. In this 
section a careful discussion of these formal power series solutions is presented. 
The analysis is generalized to vacuum spacetimes with positive cosmological 
constant in all dimensions. The expansion for perfect fluid spacetimes given in 
[T5| is also revisited. 

Consider the vacuum Einstein equations with cosmological constant A for a 
spacetime of dimension n + 1 with n > 2. An n + 1 decomposition with lapse 
equal to one and vanishing shift results in the constraint equations 



R 



k ab k ab + (trfc) 2 = 2A 
V a fc a 6 - V fc (trfc) = 



(1) 
(2) 



and the evolution equation 



d t k\ 



R a b + {tik)k a b 



2A 



(3) 



n- 1 
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Here g a b is the spatial metric with Ricci tensor R a j, and scalar curvature R, k a b 
is the second fundamental form and indices are raised and lowered using g a b 
and its inverse. Let a a b be the trace-free part of the second fundamental form 
and R a b the trace-free part of the spatial Ricci tensor and define the following 
quantities: 

E a b = d t a a b-[R a b + {trk)a a b } (4) 



E = d t (tvk) - 



, x9 2nA 
R+ (trfc) - 



(5) 



n- 1 

C = R - k ab k ab + (trfc) 2 - 2A (6) 
C a = V b k b a - Va(trfc) (7) 

The Einstein equations are equivalent to the vanishing of the evolution quantities 
E a b and E and the constraint quantities C and C a . These quantities are linked 
by the following consistency conditions 

d t C = 2(trfc)C -2V a Ca~2a a b E b a + 2(l-l/n){trk)E (8) 
d t C a = (trk)C a ~(l/2)V a C + V b E b a -(l-l/n)V a E (9) 

To do a perturbative analysis of equations CO-10 consider the formal power 
series 

IHti i 1 -Ht i 2 ~2Ht , 3 -3Ht , \ / 1 n'i 

9ab = e (g ab + g ab e + g ab e + g ab e + . . .) (10) 

where H is a constant. The n + 1 form of the Einstein equations are imposed 
on this expression in a suitable sense. It turns out that for consistency it is 
necessary to choose H = \/2A/ (n(n — 1)) and so, in particular, A must be 
positive. Products of formal series are defined in the obvious way that the 
terms in the individual series are multiplied and the resulting terms with the 
same power of e~ Ht collected. The derivatives of a formal series with respect 
to the space and time variables are defined via term by term differentiation. In 
order to impose the Einstein equations it is also necessary to have a definition of 
the inverse g ab of the formal power series metric g a b- This can be done uniquely 
by requiring that the relation g a b9 hc = 5^ holds. This allows the coefficient of 
order m in the series for g ab to be expressed in terms of the coefficients in the 
series for g a b up to order m. Setting g° ab = 5 ab and g™ b = for m > gives an 
exact solution of the Einstein equations. In the case n = 3 it is the de Sitter 
solution ( 7 1, p. 125). 

Given any tensor T, let (T) m denote the coefficient of e~ mHt in the expansion 
of T. With this notation (g a b) m = 5"fc +2 - lt f°ll° ws lrom GDI that (R a b ) m = 
for m = and m = 1. It also follows directly from H10JI that (trfe)o = —nH 
and (cr a b)o = 0. This is consistent with the vanishing of the coefficients of 
all evolution and constraint quantities for m = 0. The vanishing of (E)\ and 
{E a b)i implies that (trfc)i = and (cr%)i = 0. It follows that (C)i = and 
(Ca)i = and this ensures the consistency of the series up to order m = 1. 
Using the relation dtg a b — —2gack c b shows that g x ab — and this in turn implies 
that {R a b )3 = 0. 
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The relations between coefficients for m > 2 will now be written down. The 
summation indices p and q in the following formulae are assumed to be no less 
than two. The evolution equations © imply the recursion relations 



(n-m)(a a b ) m =H- 1 



(<J a b) P (trk) q + (R a b ) r 



j+q—m 



and 



(2n-m)(trfc) m = H^ 1 
The Hamiltonian constraint gives 



(trfc)p(trfc), + (R) 7 



p-\-q=m 



(2n-2)(trk) rn = H- 1 J £ [-(fc a 6 ) P (fc b Q ) g + (trfc) p (trfc),] + (i?) m 1 

\p+q=m ) 



(11) 



(12) 



(13) 



and the momentum constraint J3J) gives 

V° Q (fc Q b ) m -Vg(trfc) m = £ [(r= b ) p (fc a c ) g - (T a ac ) p (k c b) q ] (14) 



Here V° is the covariant derivative associated to g® b . The consistency conditions 
(|SJ| and (PJ relating evolution and constraint quantities imply that if (E)k, 
(E a b)k, (C)fe an d (C a )k all vanish for k < m — 1 then 



(2n-m)(C) f 
(n - m)ff(C ) T 



-2(n-l)(£) r 



(15) 



-(l/2)V a (C) m + V°(£ b 

a)m 

(1 - l/n)V a (i;) m (16) 



Consider first the case n — 3. The form IjlUI) of the series for g a i,, taking 
account of the vanishing of the coefficient g^ b , is contained in [18|. The following 
theorem formalizes some of the statements in ^SJ. Here smooth means C°°. 

Theorem 1 Let A ab be a smooth three-dimensional Riemannian metric and B ab 
a smooth symmetric tensor which satisfies A ab B a i, = and V°-B a b = 0, where 
the covariant derivative is that associated to A ab . Then there exists a unique 
formal power series solution of the vacuum Einstein equations with cosmological 
constant A > of the form (fTU|) with g® b = A ab and g^ b = B ab . The coefficients 
g™ b are smooth. 

Proof The coefficients (k a t,) m determine the coefficients (g a b)m recursively. For 
substituting lfTU|) into the relation d t g a b = —2g ac k c b gives 



mHgZ = 2g Q ac {¥ 



b ) m + 2 J2 9 P ac(k c b) q 

p-\-q—m 



(17) 



and hence an equation which expresses (g a b)m-2 in terms of (k a b) m and lower 
order terms for any m > 2. Thus in order to prove the theorem it is enough to 
show that equations (|ll|) - (|14fl determine the coefficients (k a b) m uniquely and 
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that when the coefficients have been fixed in this way all the equations (|f f |> — (II 411 
are satisfied. The coefficient (k a b) m is determined by (|f f \ and l|f 2fl for all to > 2 
except to = 3 and to = 6. The coefficient (k a b)3 is determined by using the 
condition that g^ b = B a b. The coefficient (<J a b)Q is determined by (|lf \ while 
(trfc)6 is determined by JT3J. By construction the evolution equation (|f 2fl is 
satisfied for all values of to except possibly to = 3 and m = 6 while l|f f t is 
satisfied except possibly for to = 3. The fact that B a b has zero trace ensures 
that (|f 21) is satisfied while l|f f t is automatic for to = 3. ft will now be shown 
by induction that ljTT|) - ((r4"Jl hold for all m. Except for m = 6 only (|T3|) and iffi)! 
need to be verified. The equations (|f f |> — (II 4|) hold for to = f. For 2 < to < 5 
the inductive step from m - 1 to m can be carried out as follows. When to ^ 3 
the consistency condition (|f 51) shows that (C) m = and then the consistency 
condition (|TB|l shows that (C a ) m = 0. That (C)s and {C a )s are zero follows 
from the conditions on B a b in the hypotheses of the theorem. Knowing that 
the equations (| 1 f |) - ljf 4|) hold for to < 5 implies using (|f 5JI that (£% = 0. By 
construction (C)6 = and then it is straightforward to obtain (C a )Q — from 
<|I 6|) . For m > 7 we can proceed as for 2 < m < 5. This completes the proof. 

Remark If P a b and P denote the Ricci tensor and Ricci scalar of g Q ab respectively, 
then g\ h = H- 2 (P ab - (l/A)Pg° ab ), a relation given in QSj. 

The theorem just proved can be generalized directly to all larger odd values 
of n, as will now be shown. 

Theorem 2 Let A a b be a smooth n-dimensional Riemannian metric with n odd 
and B a [, a smooth symmetric tensor which satisfies A ab B a b = and \7 a B a b = 0, 
where the covariant derivative is that associated to A a b- Then there exists 
a unique formal power series solution of the vacuum Einstein equations with 
cosmological constant A > of the form fTOjl with g® b = A a b and g™ 6 = B a b- 
The coefficients g™ h are smooth. 

Proof Let s be an integer such that 2s + 1 < n and (k a b) m — for all odd to 
with m < 2s — 1. If follows that g™ h = for all odd m with to < 2s — I and that 
(9 ab )m = for all odd m with to < 2s+l. Putting this information into the Ricci 
tensor shows that (R a b) m = for all odd m with m < 2s + I. Then (|II|) and l|I2|l 
imply that (k a b )2s+i — 0. It can then be proved by induction that (k a b ) m = 
vanishes for all odd to with to < n. From this point on the proof is very 
similar to that of the previous theorem. The coefficients (k a b) m are uniquely 
determined for all values of to except m = n and m = In. The coefficient (k a b)n 
is determined by using the condition that = B a b- The coefficient {<J a b)2n 
is determined by lfTT|l while (trfc)2n is determined by fl^jl . By construction the 
evolution equation l|I2|l is satisfied for all values of to except possibly to = n 
and to = 2n while Ullfl is satisfied except possibly for m = n. The fact B a b has 
zero trace ensures that l|I2() is satisfied while (|11[) is automatic for to = n. These 
statements make use of the fact that the odd order coefficients of k a b of order 
less than n vanish. It will now be shown by induction that (| 1 1 1) - (f 1 4|) hold for 
all to. Except for m = 2n only l|I3|) and (|I4|I need to be verified. The equations 
l|II|) - (|I4ll hold for to = I. For 2 < to < 2ti — 1 the inductive step from to — I 
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to m can be carried out as follows. When m ^ n the consistency condition 
<|15|) shows that (C) m = and then the consistency condition (|16fl shows that 
(Ca)m = 0. That (C)n and (C a ) n are zero follows from the conditions on B ab 
in the hypotheses of the theorem. Knowing that the equations <|11|) - (|14|) hold 
for m < 2n — 1 implies using (|15|l that (E)2 n = 0. By construction (C)2n — 
and then it is straightforward to obtain (C a )2n = 0. For m > In + 1 we can 
proceed as for 2 < m < In — 1. This completes the proof. 

The case where n is even is more complicated. The form (|1U|I of the metric 
must be generalized to 

oo L m 

9a b = e 2Ht (g° ab + ]T Y,(9a b ) m ,it l e- mHt ) (18) 

m=0 1=0 

where L m is a non-negative integer for each m and L m — for m < n. Given any 
tensor T with an expansion of the above type let (T) m j denote the coefficient 
of t l e~ mHt . As before when manipulating series they are differentiated term 
by term. The recursion relations for the expansion coefficients coming from 
the evolution equations generalize as follows, where the terms not written out 
explicitly are lower order in the sense that they can be expressed in terms of 
the coefficients of k a b with m smaller: 

(n - m)(a a b ) mJ + H'^l + l)(a a b ) m , l+1 = ... (19) 
(2n - m)(trfe) ro ,; + H'^l + l)(trfc) m , m = . . . (20) 

The terms on the right hand side not written out are obtained from the terms 
on the right hand side of equations (fTTf) and (fT*S|) if the indices m, p and q 
are replaced by the pairs {m,l), (p,h) and (q,h), summing over h + h = I- 
The recursion relations implied by the constraints are identical except for the 
addition of an extra index I. In a similar way, the consistency conditions lead 
to 

(2n-m){C) m . J +H- 1 (l + l)(C) mM1 - -2(n - l)(E) m ,i (21) 
(n-m)H(C a ) m ,i+H-\l + l)(C a ) m , l+1 = -(l/2)V a (C)„ M 

+y Q b (E b a) m ,i - (1 - l/n)V„0E) m ,j (22) 

assuming that (E)k t i- (-E a b)fc,i, {C)k,i and {C a )k,i vanish whenever k < m — 1. 

By using the above relations it is possible to express go b {g a b)n-2,o as a 
function of g^ b and its spatial derivatives. We denote this schematically by 
ffo b (ffafc)n-2,o = Z(g°). Similarly, it is possible to write V a (5 Q6 ) n _ 2 ,o = Z b (g°). 

Theorem 3 Let A ab be a smooth n-dimensional Riemannian metric and B ab a 
smooth symmetric tensor which satisfies A ab B ab = Z(A) and \7 a B ab = Z b (A), 
where the covariant derivative is that associated to A ab . Then there exists a 
unique formal series solution of the vacuum Einstein equations with cosmological 
constant A > of the form l(TH)l with g® b = A ab and (g a b)n-2,o = B ab . The 
coefficients {g a b)m.i are smooth. 
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Proof As a general principle, when determining coefficients for fixed m we start 
from / = L m and proceed to successively lower values of I. If n is odd then the 
existence follows from Theorem 2. Uniqueness for n odd in the wider class 
being considered in this theorem in comparison with Theorem 2 is obtained by 
a straightforward extension of the argument given in the proof of the latter. 
Consider now the case where n is even. By analogy with the proof of Theorem 
2 it can be shown that (k a b) m ,i — for m odd. For m < n the coefficients 
are uniquely determined. The assumption that L m = in this range is thus 
unavoidable. The coefficients (trfc) nj ; are uniquely determined by l|2(JI) and vanish 
for I > 0. The coefficients {<J a b)n,i are uniquely determined for I > 1 and vanish 
for I > 1. The choice of B ab determines {a a b)n,o- For n < m < 2n equations 
and i|2*U)l determine (k a b) m ,i- Equation i|TT)|) is used to determine (cr a b)2n,i 
while the analogue of ifHIj) is used to determine (trk) 2n ,i- For m > 2n ljT§)l and 
lf2T1|) can be used again. That all field equations are satisfied at all orders can be 
proved much as in the proof of Theorem 2, always proceeding in the direction 
of decreasing / for each fixed m. 

A question left open by Theorem 3 is whether it can ever happen that any 
of the coefficients with I > are non-zero. This is equivalent to the question 
whether (a a b) n ,i is ever non-zero. It follows from the proof of the theorem that 
this coefficient is uniquely determined by g® b . In the case n = 2 the coefficient of 
interest vanishes due to the fact that the Ricci tensor of a two-dimensional metric 
is automatically traceless. For all even dimensions greater than two there are 
choices of g® b for which (a a b) n ,i docs not vanish. In fact this is the generic case. 
The coefficient of interest can be written as a polynomial expression in H . If 
there were no logarithmic terms for a given choice of g® b then all terms in this 
polynomial would have to vanish. The coefficient of H~ n+1 , which is the most 
negative power of H occurring, is a non-zero constant times P°b(trP) fc_1 , where 
k = n/2 and P a b and trP are the tracefree part and trace of the Ricci tensor 
of g® b . There are only two ways in which this coefficient can vanish. Either 
the scalar curvature of g® b vanishes identically or g® b is an Einstein metric. A 
necessary condition for the absence of logarithmic terms has now been given 
but it is unlikely to be sufficient. The coefficients of other powers of H have to 
be taken into account in order to decide this issue. 

In |18j the expansions obtained for vacuum spacetimes were extended to the 
case of a perfect fluid with pressure proportional to energy density. Formalizing 
these considerations leads to a theorem generalizing Theorem 1 above. The 
notation here is as follows: p = T 00 , j a = T 0a and S ab = T ab . The proper 
energy density and pressure of the fluid are denoted by p and p respectively, so 
that 

T af3 = O + p)u a u f3 + pg af3 (23) 

The equation of state is taken to be p — (7 — l)p with 1 < 7 < 2. In the case 
with matter evolution and constraint quantities can be defined by 

E a b = d t a a b-lR a b + {trk)a a b-8irS a b } (24) 
E = <9f(trfc) -[R+ (trfc) 2 + 47rtrS - 12vrp - 3A] (25) 
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C = R- k ab k ab + (trfc) 2 - 167rp - 2A (26) 
C a = V b k b a -V a (trk)-87r Ja (27) 

so that their vanishing is equivalent to the Einstein equations. These satisfy 
the consistency conditions (JSJ and 10 as in the vacuum case. The components 
of the energy-momentum tensor can be expressed in terms of the fundamental 
fluid variables as follows: 

p = M(1 + 7M 2 ) (28) 
f = 7 p(l + | u | 2 )W (29) 
S\ = p[ 1U a u b + ( 7 - l)<S fc Q ] (30) 

where |it| 2 = g a bU a u b . The following relations will be useful: 

d tP -(tvk)p-^(tvk)tvS = -V a j a + a\S b a (31) 

dd a -|(trfc)i° = -V b S\ + 2a\j b (32) 

It is possible to express p and u a in terms of p and j a . To see this note first 
that \j\ = 7/i(l + |m| 2 ) 1//2 |m| and that as a consequence: 

\j\ 2 /p 2 = 1 2 (l + \u\ 2 )\u\ 2 /(l + 1 \u\ 2 ) 2 (33) 

If/(x) = 7 2 a; 2 (l + a; 2 )(l + 7.x 2 )- 2 then f(x) = 2 1 2 x(l + ~ f x 2 )- 3 {l + {2- 1 )x 2 ) > 
0. It follows that the mapping from the interval [0, oo) to the interval [0, 1) 
defined by / is invertible and \u\ 2 can be expressed as a smooth function of 
\j\ 2 /p 2 for p > 0. Since p can be expressed as a smooth function of p and \u\ 2 it 
follows that it is a smooth function of p and j a . Similarly the fact that u a can 
be expressed as a smooth function of p, \u\ 2 and j a implies that u a is a smooth 
function of p and j a . 

Next Theorem 1 will be generalized to the case with perfect fluid. The 
solution is sought as a formal series where each tensor occurring is written as a 
sum of exponentials. The exponents are taken from an increasing sequence of 
real numbers M = {rrii} which tends to infinity as i — > oo. The solution is of 
the form 



Jab 



= E (9ab) mi e- miHt (34) 

m.eM 



mteM 

E Wm,^' 
mi£M 



Let integers fci, &2, ^3 and k^ be defined as follows. For 7 < 4/3 we have 
k\ = 37, &2 = 5 — 37, ^3 = 37, ^4 = 5 while for 7 > 4/3 we have k\ = 27/(2 — 7), 
^2 = (6 — 4ry)/(2 — 7), &3 = 4 and £4 = 5. In order to organize the coefficients 
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it is useful to define the relative order rhi of a coefficient of order m, . For the 
quantities g ab , k a b , p, u a , p and j a these are defined by rhi = rrii + 2, rhi = m i, 
rhi = mi — k\, rhi = mi — k%, tth = mi — ks and rhi = ffH — k± respectively. 

Theorem 4 Let A ab be a smooth three-dimensional Riemannian metric and 
B ab a smooth symmetric tensor, po a smooth positive real-valued function and 
Uq a smooth vector field. Suppose that A ab B a b = — {8n/3H 2 )po for 7=1, 
A ab B ab = for 7 > f and V a B ab = V b {A ac B ac ) + {16TT"//3H)p A bc u% where the 
covariant derivative is that associated to A ab . If 7 > 4/3 suppose furthermore 
that u a is nowhere vanishing. Then there exists a unique formal power series 
solution of the Einstein-Euler equations with cosmological constant A > and 
equation of state p = (7 — l)p, 1 < 7 < 2, of the form l|34|) with (g a b)-2 = A ab , 
(ffafc)i = B ab , = p and (u a )k 2 = Uq. The coefficients of the series are 

smooth. They satisfy (p) mi = for m, < k%, (u a ) mi = for to,; < /c 2 and, 
except for m, = —2, the coefficient (g a b)m.i vanishes for m, < 0. 
Proof Consider a formal series solution whose coefficients vanish in the ranges 
indicated in the statement of the theorem. With the given values for k\ and ki 
the matter terms do not contribute to the equations for the coefficients of k a b 
below order three and thus all statements made about these coefficients in the 
vacuum case can be taken over without change. This follows from the fact that 
p, S a b and j a are all 0(e~ 3Ht ). The exponent in this estimate can be improved 
except in the case of p with 7 = 1 and in the case of j a with general 7. 

The proof splits into several cases. Suppose first that 1 < 7 < 1/3. Then 
|u| = o(l) and so in leading order p = p 1 j a = 7^|w|u a and S a b = (7 — 1)^5^. 
Thus the following relations are obtained: 

(m l -3j)(p) mi = ... (35) 
(mi-5)(j a ) mi = ... (36) 

The terms not written out explicitly are lower order in the sense that they are 
combinations of terms of lower relative order than m,. There is one subtlety 
involved in showing this. In the case 7 = 1 the expression V b S a b gives rise to a 
term which, looking at the exponents, is not lower order. However the coefficient 
of this term contains a factor 7 — 1 and so the term vanishes for 7 = 1. The 
Einstein equations give: 

(3-m i )(<7 a 6 ) ww = ... (37) 

and 

(6 - mi)(tvk) mi = -12ttH- 1 (p)„ h + ... (38) 

The terms on the right hand side of the last two equations not written out 
explicitly are lower order. The one explicit term on the right hand side of the 
last equation is also lower order except in the case 7 = 1. The energy- momentum 
quantities p and j a are linked to the matter quantities p and u a by the relations 

(P)m ( = Pnn + ■■■ (39) 

(u a ) mi = 1 - 1 (p- 1 f) mz + ... (40) 
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Fix a value of rhi and suppose that all coefficients with lower relative order 
have been determined. Consider the equations for (p) TOi and (u a ) mi . These 
coefficients are determined uniquely unless rhi = and if m, < they vanish. 
When rhi = they are determined by the conditions on (p)^ and (u a )k 2 in the 
hypotheses of the theorem and the equations relating p and j a to p and u a . The 
latter relations also fix the coefficients of p and u a of the given relative order 
when rhi > 0. Next consider the equations for (trfc) mi and (k a b)mi- By what has 
been said above we may assume that rrii > 3. The unique determination of the 
coefficients of the given relative order can be shown using the same procedure 
as in the vacuum case. The additional terms are either already of lower relative 
order, and hence known, or have been determined in the proceeding discussion of 
the matter equations. By induction on i it can be concluded that all coefficients 
are uniquely determined. The fact that all field equations are satisfied can 
be shown much as in the vacuum case since the compatibility conditions are 
identical. 

Now consider the case 4/3 < 7 < 2. The assumption that Uq is nowhere 
vanishing implies in this case that = o(l) and in leading order p = jp\u\ 2 , 

j a = -fp\u\u a and S a b = -fpu a Ub- The following relations are obtained: 

(rm - 4)(p) mi = ... (41) 
K-5)(.f) roi = ... (42) 

The Einstein equations give the same relations as in the previous case. For 
4/3 < 7 < 2 the energy-momentum quantities p and j a are linked to the matter 
quantities p and u a by the relations 

Mm* = (2-7)- 1 (p(l-blVp 2 ))m i +--- (43) 

(u a ) mi = ((2- 7 )/7) 1/2 (p- 1 (l-b1 2 /p 2 )" 1/ V) mi +... (44) 

Using these facts we can proceed as in the case 1 < 7 < 4/3. 

Consider finally the case 7 = 4/3 where \u\ tends to a finite limit, in general 
non-zero, as t — > 00. The difference in comparison to the cases already treated 
is that the relations between p and j a on the one hand and p and u a on the 
other hand cannot be inverted explicitly in leading order. However the fact, 
shown above, that the relevant mappings are invertible has an equivalent on the 
level of formal power series. For if / is a smooth function between open subsets 
of Euclidean spaces then f(x + y) can be written formally in terms of a Taylor 
series about x. The resulting expression contains the derivatives of / evaluated 
at x multiplied by powers of y. If y is replaced by a formal power series without 
constant term then a well-defined formal power series for f{x + y) is obtained. 
Thus the same method can be applied as in the previous cases, allowing the 
proof of the theorem to be completed. 

A case which has been excluded in the above theorem is that where 7 > 
4/3 and u a may vanish somewhere. In that case the kind of series which has 
been assumed in the theorem is not consistent. For if it is assumed that an 
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expansion of this kind is possible this leads to different rates of decay for certain 
quantities, for instance fi, depending on whether |u| does or does not vanish. 
As a consequence V a /x//z will be unbounded as t tends to infinity although /i 
is nowhere zero. This contradicts the assumptions which have been made. The 
situation is reminiscent of the spikes observed near the initial singularity in 
Gowdy spacetimes ^7] and so we may speculate that in reality inhomogeneous 
features develop in \i so that the density contrast blows up as t — > oo. This 
behaviour for 7 > 4/3 is not consistent with the usual picture in inflationary 
models where the density contrast remains bounded at late times. The issue 
deserves to be investigated further. 

It is interesting to ask whether the expansions for a fluid presented here 
can be extended to the case of collisionless matter. If they can then the result 
probably resembles that for dust. Limited expansions in some special cases are 
already known 

Note that the analysis of vacuum spacetimes in this section has a close ana- 
logue for Riemannian (i.e. positive definite) metrics. A solution of the Einstein 
equations with positive cosmological constant in the Lorentzian case corresponds 
to an Einstein metric with negative Einstein constant in the Riemannian case. 
The equations obtained for a positive definite metric are 

- R - k ab k ab + (trfc) 2 = -(n-l)K (45) 
V a k a b - V 6 (tr*) = (46) 
8 t k a b = -R a b + (trk)k a b + 2KS a b (47) 

where K is the Einstein constant, i.e. the n + 1-dimensional metric satisfies 
R a p = Kg a p. Asymptotic expansions for this case have been investigated in 
the literature on Riemannian geometry |3] and string theory 



3 From minimal to full asymptotics 

In the last section consistent formal asymptotic expansions were exhibited for 
a number of problems. In this section it is shown that minimal information 
about the asymptotics implies the full expansions given in the last section. For 
simplicity we restrict consideration to the vacuum case. The following lemma 
will be used: 

Lemma 1 Consider an equation of the form 

5 tW + fc U = ^ Um / e - mt + 0( e - Jt ) (48) 

in. I 

for a vector- valued function u(t), where j ^ k and m < j in the sum. Then 
there are coefficients u m j, m < j, such that 

u = J2^mjt l e- mt + 0(e-^) (49) 

rnj 
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If 148f) may be differentiated term by term with respect to t as often as desired 
the same is true of (|49ll . 

Proof Note first that dt(e kt u) is equal to a sum of explicit terms with a term 
of order e' fc_: '- ) *. Each of the explicit terms has an explicit primitive which is a 
sum of terms of the same general form and the same value of m but in general 
several values of I . Thus we can absorb these terms into the time derivative and 
write 

d t (e kt (u - u mtl t l e- mt )) = 0(6^-^) (50) 

m.l 

with m < j in the sum. If j < k we can integrate this relation directly to get the 
desired result. If j > k then the expression which is differentiated with respect 
to time converges to a limit as t — > oo, which can be called Uk,o- This gives the 
desired result in the latter case. 

If the assumption on time derivatives is satisfied then dtu satisfies an equa- 
tion of the same form as that satisfied by u. Hence dtu has an asymptotic 
expansion 

d t u = ™ m ,lt l e- mt + 0(e-^) (51) 

rn .1 

Integrating this from to to t and using Q49J) gives 

[ ™m,is l e~ ms ds = C + J2 u m ,lt l e- mt + 0(e-^) (52) 

m,l *° m,l 

for a constant C. It follows that the coefficients w m ,i are obtained from u m> i 
by term by term differentiation. This process can be repeated for higher order 
derivatives with respect to t. 

Remark If the quantities in (|48|) depend smoothly on a parameter and the 
equation may be differentiated term by term with respect to the parameter 
then the same is true for the solution. 

Theorem 5 Let a solution of the vacuum Einstein equations with cosmological 
constant A > in n + 1 dimensions be given in Gauss coordinates. Suppose 
that e~ 2Ht g a b, e 2Ht g ab , e 2Ht a a b and their spatial derivatives of all orders are 
bounded. Then the solution has an asymptotic expansion of the form given in 
Theorem 3. The expansion remains valid when differentiated term by term to 
any order. 

Proof The Hamiltonian constraint can be used to express trfe in terms of the 
scalar curvature R, a a b and A, giving 



trfe 



n- 1 



[-R + <J a b <j b a ) + n 2 H' 



1/2 



(53) 



It follows from the assumptions of the theorem that trfc = — nH + 0(e 2Ht ) 
and that this relation may be differentiated term by term with respect to the 
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spatial variables. Now 

d t (e~ 2Ht g ab ) = -2e- 2Ht g ac {k\ + HS c b ) (54) 

The right hand side of this expression is 0(e~ 2Ht ) and so there is some g® b such 
that 

e- 2Ht g ab = g ab + O(e- 2Ht ) (55) 

and corresponding relations hold for spatial derivatives of all orders. Using the 
evolution equations it can be seen that these relations can also be differenti- 
ated repeatedly with respect to time. The proof now precedes by induction. 
The inductive hypothesis is as follows. There exist coefficients (g a b)m,i and 
(fc a h)m,/,0 <m< M such that 

M L m 

g a b = e 2fft (^ *^2(gab)m-2,it l e~ mHt + (jab) = [g a b]M + e 2Ht g a b (56) 
m=0 i=a 
M L m 

m=0 1=0 

where g ab and k a are 0(e~( M+ ^ Ht ) and similar asymptotic expansions hold for 
all derivatives of these quantities. Here e is a constant belonging to the interval 
(0, 1). The inductive hypothesis is satisfied for M = 1. If these expressions are 
substituted into the Einstein equations then the expansion coefficients written 
explicitly satisfy the same relations as in the analysis of formal power series 
solutions carried out above. It is convenient to write the evolution equations in 
the following form: 

d t g a b = -2g ac a c b - (2/n)(trfc + nH)g ab (58) 
d t a a b + nHa a b = (trfc + nH )a a b + R a b (59) 
dt(txk + nH)+2nH(txk + nH) = (trfc + nH) 2 + R (60) 

where g ab — e~ 2Ht g ab . Using the inductive hypothesis it follows that if each 
quantity Q in these equations is replaced by the corresponding quantity [Q]m+i 
then equality holds up to a remainder of order e -( M + 1 + i ) Ht j n (|59|| and (|60|) . 
Using this information shows that the corresponding statement holds in l|58|l 
with a remainder of order e -\ M - l +v Ht _ Thus the quantities [Q]m+i satisfy a 
system of the type occurring in Lemma 1. It follows from that lemma that the 
inductive hypothesis is satisfied with M replaced by M + 1. 

In Jl] results similar to those of this section were obtained using a different 
coordinate system. The time coordinate used there satisfies the condition that 
that the lapse function is proportional to the inverse of the mean curvature of 
its level surfaces. This means that the foliation of level surfaces is a solution 
of the inverse mean curvature flow, a fact which raises serious doubts whether 
such coordinates exist in forever expanding cosmological spacetimes, as will 
now be explained. The inverse mean curvature flow for hypersurfaces is defined 



13 



by the condition that a hypersurface flows with a speed equal to the inverse 
of its mean curvature in the normal direction. In the case of a Riemannian 
manifold it was used in the work of Huisken and Ilmanen ^U] on the Penrose 
inequality. For spacelike hypersurfaces in a Lorentzian manifold it was studied 
in [5]. If a spacelike hypersurface with positive expansion (i.e., in the convention 
used here, with trfc < 0) is given then there is a local solution of the inverse 
mean curvature flow in the contracting direction. Moreover, under reasonable 
assumptions on the nature of the singularity, there is a global solution. In the 
expanding direction, in contrast, the equation is backward parabolic and it is 
to be expected that there is no local solution for general initial data, i.e. for 
a general starting hypersurface. This is an analogue of the fact that the heat 
equation cannot be solved backwards in time. 

4 Relations to conformal infinity 

There is a relation between the expansions discussed in the last two sections 
and the concept of conformal infinity. In this section only the Einstein vac- 
uum equations are considered. Define T = H~ 1 e~ Ht . Then spacetime metric 
corresponding to (fT%)) becomes 

(HT)~ 2 [—dT 2 + (g° ab + J2(9 ab )m,l(-l) l H- l (log(HT)) l (HTr)} (61) 

m=0 1=0 

It is conformal to a metric which is non-degenerate at T — and is written 
in Gauss coordinates. If there are no non-vanishing coefficients with / > the 
conformal metric (or unphysical metric) is smooth at T — 0. This is for instance 
the case when n is odd. 

In the case n = 3 Friedrich 0], |S] has used conformal techniques to prove 
results which, as shown in the following, imply that spacetimes evolving from 
initial data close to standard initial data for de Sitter space indeed have asymp- 
totic expansions of the type presented in the last section. The method used, 
based on the conformal method, is only known to work in the case n = 3. The 
occurence of logarithms in the expansions for even values of n cast doubt on 
the possibility of implementing an analogous procedure in that case. There are 
also problems for n — 3 if matter is present. For conformally invariant matter 
fields the method can be used but for other types of matter, e.g. a perfect fluid 
with linear equation of state, there is no straightforward way of doing this. The 
non-integer powers occurring in the formal expansions for this case make the 
application of the method problematic. Note, however, that a similar problem 
has been overcome in the study of isotropic singularities [2]. 

Consider the de Sitter solution with a slicing by intrinsically flat hypersur- 
faces, as described in Section 2, with the slicing being given by the hypersurfaces 
of constant t. We may assume for convenience that the solution has been iden- 
tified in a way which is periodic in the spatial coordinates. Consider initial data 
which is a small perturbation of the initial data induced by this model solution 
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on a hypersurface of constant time. The smallness can be measured in the sense 
of uniform convergence of a function and its derivatives of all orders. Then, 
according to section 9 of |5] , the perturbed solution has a Cauchy development 
which is asymptotically simple in the future. This means that the solution g a p 
is conformal to a metric g a p = fl 2 g a p with fl > in such a way that g a p and 
fl have smooth extensions through a hypersurface where fl vanishes. We may 
choose coordinates in the unphysical metric in the following way. Set T = fi 
and choose the spatial coordinates to be constant along the curves orthogonal 
to the hypersurfaces of constant f2. In these coordinates the conformal metric 
takes the form 

- H- 2 a 2 df 2 + g ab dX a dX b (62) 

where a is a function of T and X a . The condition — 3V a f2V a £l = A (see Lemma 
9.2 of [5]) implies that a = 1 for T = 0. In order to compare this with the 
expansions in the previous sections we need to transform to Gauss coordinates 
with respect to the physical metric g a p. As a first step let T — e~ HT . Then the 
physical metric becomes 

- a 2 dT 2 + g ab dX a dX b (63) 

with g a i, = e 2HT g ab . The following lemma shows that Gauss coordinates of 
a suitable kind can be introduced. The hypotheses make use of the following 
inequalities 

I Sab | 

la^l + IIXJ 
|a-l| + |0ra| + |0„a| 
\k a b \ + \ti-k + 3H\ + \g ab \ 

The metric l|63|) above satisfies inequalities of this type together with corre- 
sponding inequalities for spatial derivatives of all orders. In fact the estimates 
for k a b and trfc + H are only obviously satisfied with the bound Ce~ HT . How- 
ever this can be improved by using equations (|59|) and l|6U|) at the end of the 
last section, or rather their equivalents in the presence of a non-trivial lapse 
function. 



< Ce 2HT (64) 

< C (65) 

< Ce- HT (66) 

< Ce- 2HT (67) 



Lemma 2 Consider a metric of the form (|63|l on a time interval [To, oo) and 
assume that there is a constant C > such that the inequalities l|64l) - (|67|l are 
satisfied, together with the corresponding inequalities for spatial derivatives of 
all orders. Then for To sufficiently large there exists a Gaussian coordinate 
system based on the hypersurface T = Tq which is global in the future. The 
transformed metric satisfies the hypotheses of Theorem 5. 

Proof To construct Gaussian coordinates it is necessary to analyse the equations 
of timelike geodesies. In 3+1 form these are 

d 2 T / dT\ 2 „ dX a dT dX a dX b 

—r + a o T a[ — ) +2a V a a——— + a k ab — — = 

dr V dT I dT dr dT dT 
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(69) 



It is helpful for the following analysis to rewrite one of the terms: 

2 , l dX a dT nTT dX a nTT dX a ( dT , 
--a(tik) — — = 2H— — + 2H— — 1 

3 y ' dr dr dr dT \dr 

2, , nTT ^dX a dT 2, ^dX a dT^ n , 
- - (trfc + 3tf — — - - (a - 1 trfc -j- -p 70 
3 dr dT 6 dT dr 

These equations are to be solved for functions T(r, a; b ) and X a (r, x 6 ) with initial 
values T = T , dT/dr = 1, X a = x a and dX a /dr = at r = T . Strictly 
speaking Gaussian coordinates based on T = T would differ from this by a time 
translation by T but it is convenient here to work with this slight modification. 
Consider now a solution of these equations on an interval [To, r*) and suppose for 
later convenience that To > 0. There is a t* > To for which a solution does exist. 
We assume that on this interval \dX a /dT\ < Ce~ 2ffT and that \dT/dT - 1| < e 
for some e G (0, 1/3). For given C and e there exists an interval of this kind. 
On this interval e~ T < e~ eTo e~( 1 ~ £ ' T , e T ~ T < e £ ^ T_T °' and inequalities of the 
following form hold, where C is a positive constant depending only on C and 
e. 

^ = f(r), \f(T)\<Ce-WT (71) 

^T + 2H*£- = g(T), \g(T)\<C'e-^ (72) 

It follows from the first of these that 

\dT/dT-l\ <C'e- (1 -^ HTo (73) 

For To large enough this strictly improves on the estimate originally assumed for 
dT/dr—1. For small e the quantities dX a /dT can be seen to decay exponentially 
with an exponent which is as close as desired to —2. The fact that we are dealing 
with timelike geodesies parametrized by proper time leads to the relation 

2 (dT\ 2 dX a dX b 

This implies that \dT/dr — 1| = 0(e~ HT ). Putting this back into the evolution 
equation for dX a /dr shows that it is 0(e~ 2HT ). By choosing To large enough 
the decay estimate for this quantity is recovered and in fact strengthened. Con- 
sideration of the longest time interval on which the original inequality holds 
shows that t* — oo. The estimates we have derived up to now hold globally. 
The estimate for dT/dr obtained above implies that there are positive constants 
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C\ and Ci such that C\t < T < C2T. Hence in estimates we can replace e~ T 
by e~ r if desired. 

Next we would like to obtain corresponding estimates for the spatial deriva- 
tives of dT/dr and dX a /dr of all orders. Consider the result of differentiating 
the geodesic equations with respect to the spatial variables. This leads to esti- 
mates of the form 

This allows us to show that the first order spatial derivatives of the key quantities 
satisfy the estimates analogous to those satisfied by the quantities themselves. 
The same argument can be applied to estimate spatial derivatives of any order 
inductively. Now all the desired information about existence and decay of T 
and X a has been obtained. It remains to show that they form a coordinate 
system. This follows from the fact that the initial values of dT/dr, dX a /dr 
and dX a /dx b are one, zero and 6% respectively and the exponential decay of 
their time derivatives which has already been proved. This completes the proof 
of the lemma. 

Combining Lemma 2 and Theorem 5 shows that the spacetimes constructed 
by Friedrich admit global Gaussian coordinates in which they have an asymp- 
totic expansion of the form of (|10|1 . Hence any initial data close to that for 
de Sitter on a flat hypersurface evolves into a solution having an asymptotic 
expansion of the form given by Starobinsky. 

5 The spatially homogeneous case 

This section is concerned with spatially homogeneous solutions of the vacuum 
Einstein equations in n + 1 dimensions with positive cosmological constant. It 
is assumed that the spatial homogeneity is defined by a Lie group G, supposed 
simply connected, which acts simply transitively. We restrict to spacetimes 
such that all left invariant Riemannian metrics on G have non-positive scalar 
curvature. In the case n = 3 this corresponds to Bianchi types I to VIII. 
Information on the case n = 4 can be found in [Sj. It will be shown that the 
spacetimes of the type just specified have asymptotic expansions with all the 
properties of the formal expansions in Theorem 3. 

A spatially homogeneous spacetime of the type being considered can be 
written in the form 

-dt 2 + gij (t)e l (g>e j (77) 

where {e 1 } is a left invariant frame on the Lie group G. Basic information about 
the asymptotics of these spacetimes in 3+1 dimensions are given by Wald's the- 
orem |21| . which provides information on the behaviour of the second funda- 
mental form as t — > 00. This can easily be generalized to the present situation. 
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Using the condition on the sign of R, the Hamiltonian constraint implies that 
on any interval where a solution exists (trfc) 2 > = (nH) 2 . Combining the 
Hamiltonian constraint with the evolution equation for trfc gives 

9 t (trfc) > i(trfc) 2 - -^-A (78) 

In particular trfc is non-decreasing. These facts together show that trfc is 
bounded. Now it will be shown that trfc — > —nH as t — > oo. For 

d t {trk + nH) > -(-trfc + nH) (-trfc - nH) (79) 
> -2iT(trfc + nH) (80) 

It follows that trfc = —nH + 0{eT 2Ht ). Using the Hamiltonian constraint then 
gives a 1 i — 0(e~ 2Ht ). This bound can be used to get information on ov,- as 
in • Then it is possible to proceed exactly as in the proof of Proposition 2 
in J2j to show that e~ 2Ht gij, e 2Ht g % i and e Ht a l j are bounded. Then equation 
(58) can be used as in the previous section to improve the last statement to 
the boundedness of e 2Ht a' l j. The fact that gy, g^ and fey are bounded on any 
finite time interval implies that the solution exists globally in time. 

We are now in a situation very similar to that of Theorem 5. However the 
estimates we have are expressed in term of frame components. Choosing a 
coordinate system on some subset of the Lie group G with compact closure will 
give us uniform asymptotic expansions for the components in that coordinate 
system. Conversely uniform asymptotic expansions for the components in a 
coordinate system of this type give corresponding asymptotic expansions for 
the frame components. In this case we will say that the asymptotic expansions 
are locally uniform. If an expansion of this type holds for a given quantity it 
also holds for all spatial derivatives in the coordinate representation. 

The proof of Theorem 5 uses only arguments which are pointwise in space 
and so it generalizes immediately to the case of locally uniform asymptotic ex- 
pansions. It can be concluded that for a spatially homogeneous spacetime of 
the type under consideration locally uniform asymptotic expansions of gener- 
alized Starobinsky type are obtained. Restricting to a coordinate domain with 
compact closure uniform asymptotic expansions are obtained. In general these 
expansions will contain logarithmic terms. Consider for instance the case of the 
Lie group H x R where H is a three-dimensional Lie group of Bianchi type other 
than IX. Let the spacetime be such that the spatial metric at each time is the 
product of a metric on H with one on R. This is consistent with the constraint 
equations. For instance the initial data can be chosen to be the product of data 
on H with trivial data on R. The data are invariant under reflection in R and 
this property is inherited by the solutions. Suppose that H admits no metric 
of vanishing scalar curvature. This is the case for every Bianchi type except I 
and VIIo . Then the metric g® b corresponding to this solution has non- vanishing 
scalar curvature and is not an Einstein metric. Hence logarithmic terms are 
unavoidable. 
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6 Fuchsian analysis 



Fuchsian systems are a class of singular equations which can be used to prove 
the existence of solutions of certain partial differential equations with given 
asymptotics ^T], ^2]> |lfc>| - It will be shown that Fuchsian methods allow the 
construction of solutions of the vacuum Einstein equations with positive cosmo- 
logical constant in any number of dimensions which have asymptotic expansions 
of the type given in Section 2 and depend on the same number of free functions 
as the general solution. 

Before coming to the specific problem of interest here some general facts 
about Fuchsian equations will be recalled. The form of the equations is 

tdu/dt + N(x)u = tf(t,x,u,u x ) (81) 

Here x denotes the spatial coordinates collectively and u x the spatial derivatives 
of the unknown u(t,x). The matrix N and the function / are required to 
satisfy certain regularity conditions and N is required to satisfy a positivity 
condition. There are forms of the regularity condition adapted to smooth and 
to analytic functions. The version adapted to analytic fucntions will be used in 
the following since it is the one where the most powerful theorems are available. 
For the precise definition of regularity see 1 , where a corresponding definition 
of regularity of solutions is also given. Roughly speaking, regularity means that 
the functions concerned are continuous in t and analytic in x and vanish in a 
suitable way as t — > 0. Consider now an ansatz of the form 

oc L m 

u(t,x)=J2J2 u m,it m (logt) 1 (82) 

m=0 1=0 

By analogy with what was done in Section 2 we can ask whether the equa- 
tion (|81|l has a formal series solution of this kind. Suppose that this is the 
case. Fix M > 0. Then there exist coefficients u m j such that u — u — 
Ya? Lo™ u rn .it m (\ogt) 1 satisfies 

tdu/dt + N(x)u = t M+e f M (u) (83) 

for a regular function /m and a constant e > together with the corresponding 
relations obtained by differentiating term by term with respect to the spatial 
coordinates any number of times. Let v — t~ M u. Then 

tdv/dt + (N(x) + MI)v = t e g(t, x, v, v x ) (84) 

for a regular function g. Introducing t f ~ as a new time variable, we obtain an 
equation of the form (|81|l . If we assume that N(x) is bounded then by choosing 
M large enough it can be ensured that the matrix N(x)+MI is positive definite. 
Assuming that / and N are regular in the analytic sense the existence theorem 
of PP implies the existence of a unique regular solution v vanishing at t = 0. 
Expressing u in terms of v gives a solution of the original equation which has 
the given asymptotic expansion up to order M. 
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Consider now the slightly more general equation 

tdu/dt + N(x)u = tf(t, x, u, u x ) + h{u) (85) 

In order to have a consistent formal power series solution suppose that for some 
functions wo,o & n d iti^o we have Nuq.o = /i(ito.o) and 

(N + I - Dh(u 0fi ))u lfi = /(0, u 0t0 ) (86) 

Here Dh denotes the derivative of h as a map between Euclidean spaces. Sup- 
pose further that the equation admits a formal power series solution with coef- 
ficients wo,o and ui,o and Lq = L\ = 0. If these conditions hold then u satisfies 
the original equation if v = u — uq^ — Uifit satisfies a Fuchsian system and 
vanishes at the origin. Thus an existence theorem is obtained. 

To make contact with the Einstein equations we start with the equations 
<|58|) - (|t)U[) and set r = e~ Ht . Then an equation of the form l|85|) is obtained, 
with u = (g a b, & a b, trfc + nH). If it is assumed that the variables a a b and 
(trfe + nH) vanish at r = then the consistency conditions on ito,o and Mi,o 
arc satisfied. The fact that consistent formal expansions were shown to exist 
in Section 2 allows the above procedure to be carried through. If the data A a \, 
and B a b are chosen to be analytic then this gives an existence theorem for the 
Einstein evolution equations with A a b and B a b prescribed as in Theorem 3. In 
this context it is important to note that if A a b and B a b are analytic all the 
coefficients in the formal expansions whose existence is asserted in Theorem 3 
are also analytic. In order to see that a solution of the Einstein equations is 
obtained it suffices to show that the constraint equations are satisfied. Note 
that it follows from the results of Section 2 that the constraint quantities vanish 
to all orders at r = but since the solution is not analytic at r = this 
does not suffice to conclude that the constraint quantities vanish everywhere. 
To see that they do we need to write the consistency conditions (JSJ and © in 
Fuchsian form. Using the fact that the Einstein evolution equations are satisfied, 
and introducing C a = e~ fH C a , these equations can be written as 

d t C + 2HC = 2(trfc + H )C - 2e m V a C a (87) 
d t C a + 2HC a = (trk + H)C a -(l/2)e- m V a C (88) 

Setting t = e~ Ht gives a system of the form (|8*T1l and since C and C a tend to zero 
as t — ► both of these quanitities vanish as a consequence of the uniqueness 
theorem for Fuchsian systems and the constraints are satisfied. The solution 
of the Einstein equations has an asymptotic expansion of the form given in 
Theorem 3 truncated at any given finite order. Applying Theorem 5 shows that 
this solution has an asymptotic expansion of this form to all orders. The results 
obtained can be summed up as follows: 

Theorem 6 Let A a b be an analytic n-dimensional Riemannian metric and B a b 
an analytic symmetric tensor which satisfies A ab B a b — Z{A) and V a B a b = 
Zb(A), where the covariant derivative is that associated to A. Then there exists 
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an analytic solution of the vacuum Einstein equations with an asymptotic ex- 
pansion of the form (|18() with g® b = A a b and (<? Q &)n-2,o = B a b. The expansion 
may be differentiated term by term with respect to the spatial variables as often 
as desired. 

7 The wave equation on de Sitter spacetime 

In Section 4 it was shown that initial data for the vacuum Einstein equations in 
3+1 dimensions close to that for de Sitter space evolve to give a spacetime with 
asymptotics of Starobinsky type. It has not yet proved possible to obtain the 
analogous statement in higher dimensions. What is missing are suitable energy 
estimates. In this section it will be shown how a simpler model problem can 
be treated. This is the case of the wave equation Vo,V a ^> = on (the higher 
dimensional analogue of) de Sitter space. The spacetime metric in this case is 

ds 2 = -dt 2 + e 2Ht {{dx l f + ... + (dx n ) 2 ) (89) 

Written out explicitly in coordinates the wave equation takes the form: 

d 2 <j> + nHd t <f> = e- 2Ht A<j> (90) 

where A is the Laplacian of the flat metric. 

Consider the ansatz for formal solutions of the equations 

oo 

(A m (x)e~ mHt + B m {x)te- mHt ) (91) 

m=0 

Substituting this into the equation and comparing coefficients gives 

m(m - n)H 2 A m - (2m - n)HB m = AA m - 2 (92) 
m(m-n)H 2 B m = AB m _ 2 (93) 

For any n it is true that Bq = A\ = B\ = 0. In the case that n is odd assume 
that the coefficients B m vanish. Then AA m _2 = H 2 m(m — n)A m for all m > 2. 
Then it follows from Ay = that A 2 ^+i = for all integers k with 2k + 1 < n. 
The coefficients A 2 k+i with 2k + 1 > n are determined by A n . The coefficients 
A 2 k are determined by Aq. There are no further relations to be satisfied and 
so Aq and A n parametrize the general solution. If the coefficients B m are not 
assumed to be zero it can be shown that they must vanish for n odd. For n even 
we have A 2 k+i — £>2fc+i = for every positive integer k. If 2k < n then B 2 k = 0. 
For 2k < n the coefficients A 2 k are determined successively by Aq. Also B n is 
determined by Aq. Then A 2 k and B 2 k are determined for all k with 2k > n in 
terms of the coefficients already determined. Thus the general solution can be 
parametrized by Aq and A n , just as in the case n odd. The difference is that 
the series obtained contains terms which are multiples of te~ mHt . 

Let E — J (d t 4>) 2 + e~ 2Ht \V<j)\ 2 . Differentiating with respect to t and inte- 
grating by parts gives the relation dE/dt = —2HE. We can differentiate the 



21 



equation through with respect to a spatial coordinate and repeat the argument. 
This shows that all Sobolev norms of e Ht d t (f> and V<j) are bounded. By the 
Sobolev embedding theorem they and all their spatial derivatives satisfy corre- 
sponding pointwisc bounds. Thus the spatial derivatives of <fr are bounded while 
its time derivative decays like e~ Ht . As a consequence <j>(t, x) = (po(x)+0(e~ Ht ) 
for some function (j> . Comparing with the formal solutions already obtained 
we see that these estimates are not likely to be sharp. The equation for (f> is 
equivalent to the system 



Starting with the basic information on the asymptotic behaviour of <f> we already 
have the method of proof of Theorem 5 can be applied to this system. The result 
is that any solution has an asymptotic expansion of the type derived on a formal 
level above. 
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